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EFFECTS OF MAGNETIC TURBULENCE ON THE DYNAMICS OF 
PICKUP IONS IN THE lONOSHEATH OF MARS 

H. Aceves,^ M. Reyes-Ruiz & C. E. Chavez 
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RESUMEN 

Se estudian los efectos que fluctuaciones magneticas turbulentas tienen sobre 
la dinamica de iones pickup de O"*" en las regiones polares magneticas de la iono- 
funda de Martc. En particular cstudiamos su cfecto en los perfiles de velocidad 
proniedio como funcion de la altitud sobre los polos, con el fin de comparar con 
datos reciente del Mars Express; que indican que su velocidad proniedio es pequefia 
y esencialmente en la direccion antisolar. Se encuentra que, aunque las fluctuaciones 
magneticas resultan en desviaciones del simple movimiento de arrastre girotropico, 
aun amplitudes de fluctuacion mas grandes que las medidas, no son capaces dc re- 
producir los perfiles verticales de velocidad de los iones de . Conclm'mos que otro 
mecanismo fisico, diferente a uno puramente de dinamica dc particulas cargadas, 
esta actuando sobre estos iones en el terminador Marciano. Una posibilidad es que 
exista una interaccion tipo-viscosa entre el viento solar y la ionosfera Marciana a 
bajas altitudes. 

ABSTRACT 

We study some of the effects that magnetic turbulent fluctuations have on the 
dynamics of pickup ions in the magnetic polar regions of the Mars ionosheath. 
In particular we study their effect on the bulk velocity profiles of ions as a function 
of altitude over the magnetic poles, in order to compare them with recent Mars 
Express data; that indicate that their average velocity is very low and essentially in 
the anti-sunward direction. We find that, while magnetic field fiuctuations do give 
rise to deviations from simple E x B-drift gyromotion, even fluctuation amplitudes 
much greater than those of in situ measurements are not able to reproduce the 
vertical velocity profile of 0+ ions. We conclude that other physical mechanisms, 
different from a pure charged particle dynamics, are acting on pickup ions at the 
Martian terminator. A possibility being a viscous-like interaction between the Solar 
Wind and the Martian ionosphere at low altitudes. 

Key Words: Physical Processes: turbulence, magnetic fields, acceleration of parti- 
cles — Solar System: planets: Mars — Methods: numerical 



1. INTRODUCTION 

The interaction of the solar wind (SW) with So- 
lar System objects with a negligie intrinsic magnetic 
field and an atmosphere (e.g. Venus and Mars) is cur- 
rently being investigated by several space missions 
(e.g. Venus and Mars Express), that will enhance 
our understanding of plasma interaction processes in 
such environments (e.g. Ma et al. 2008). In partic- 
ular, a comparison of different models for the global 
interaction of the solar wind (SW), a weakly colli- 
sional plasma (e.g. Marsh 1994, Echim et al. 2011), 
with Mars has been done by Brain et al. (2010) 
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that included MHD, multi-fluid and hybrid models. 
All these models, as well as the gasdynamic con- 
vected magnetic field models (e.g. Belotserkovskii 
et al. 1987, Spreiter & Stahara 1994), have greatly 
advanced our knowledge of the complex interaction 
of the SW plasma and the Martian ionosphere. 

An important component of the Martian plasma 
environment, specially at low and mid altitudes, is 
the population of oxygen ions (O"*") produced by the 
UV solar radiation and charge exchange reactions. 
These are picked up by the SW convective electric 
field and removed from the planet, in turn affecting 
the solar wind interaction with Mars (e.g. Ledv- 
ina et al. 2008, Withers 2009) 
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Recently, Perez-de-Tejada et al. (2009; PdT09) 
have pointed out that the analysis of Mars-Express 
measurements, conducted with the ASPERA-3 in- 
strument (Barabash et al. 2006), indicates that the 
bulk velocity of ions, presumably of ionospheric 
origin, as they stream at low altitude over the mag- 
netic poles of the planet is mostly in the anti- 
sunward direction (see Figure[T]). This result is unex- 
pected in the context of a SW-ionosphere interaction 
in which the acceleration of pickup ions is essentially 
due to the action of the convective electric field, E. 
As shown in Reyes- Ruiz et al. [2010b (RAPIO)], in 
the context of a simplified model for the geometry 
of the IMF and SW fiow in the ionosheath of Mars, 
the bulk velocity of pickup ions in gyromotion, 
as they pass over the magnetic poles of the planet, is 
strongly dominated by the vertical component (w^, 
normal to the plane of the ecliptic) . 

A qualitatively similar conclusion to the above 
can also be reached from the results of MHD sim- 
ulations using more complicated IMF and Martian 
magnetic field models to study the dynamics of 
heavy ions in the region (e.g. Fang et al. 2008), 
and seen in earlier pickup ion trajectories in 
a draped magnetic field model (e.g. Luhmann & 
Schwingenschuh 1990, Luhmann 1990, Lichteneg- 
ger et al. 2000). In RAPIO it is argued that a largely 
horizontal bulk velocity of ions can be explained 
in terms of a viscous-likc interaction between the SW 
and ionospheric plasmas, as found in the numerical 
simulations of Reyes- Ruiz et al. (2009). 

It has been suggested that magnetic field fluc- 
tuations, as those known to be widespread in the 
ionosheath of Mars and other non-magnetic bodies 
(Tsurutani et al. 1995, Nagy et al. 2004, Grebowsky 
et al. 2004, Espley et al. 2004, Voros et al. 2008), 
may be responsible for modifying the gyromotion of 
pickup ions, enhancing pitch-angle scattering and ion 
heating and acceleration (e.g. Wang et al. 2006 and 
references therein). Najib et al. (2011) have noted 
more recently that the presence of wave activity and 
turbulence lead to wave-particle interactions that re- 
semble collisions, hence probably modifiyng ion dis- 
tribution function. In this paper we study the possi- 
bility that such modifications of simple gyromotion 
due to turbulcncc-like fiuctuations of the magnetic 
field, in the ionosheath over the magnetic poles of 
Mars, lead to a bulk velocity profile for 0+ ions as 
that reported by PdT09 from the ASPERA-3 mea- 
surements. 

It is worth emphasizing that global MHD or 
Hybrid simulations analyzing the motion of iono- 
spheric ions (Nagy et al. 2004, Fang et al. 2008, 
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Fig. 1. ASPERA-3 in situ measurements of the velocities 
profiles of SW protons {top) and pickup ions (bottom) 
over the terminator of Mars; coordinates are in a MSE 
system, and the data is taken from PdT09. Four alti- 
tudes of the spacecraft are indicated as dotted vertical 
lines. The dominance of the a;-component of the velocity 
of ions over its z component is very clear, contrary 
to what is expected from a pure U x B-drift. 

Brain et al. 2010, Kallio et al. 2010, Najib et al. 2011, 
and the review by Nagy et al. 2004,) do not include 
the effect of small-scale turbulent fluctuations of the 
magnetic held, known to exist in the ionosheath of 
nonmagnetic Solar System bodies. These models es- 
sentially use the B field obtained directly from MHD 
simulations to follow the dynamics of the 0+ ions. 
Hence, our approach is not physically modeling the 
same phenomena that the above models are able to 
investigate. On other hand, do not consider the ef- 
fect the crustal field (e.g. Acufia et al. 1998, Zhang 
& Li 2009) can have on the ion dynamics. 

This paper is organized as follows. In §2 we 
present the basic equations and methodology used 
for our analysis, including the prescription for con- 
structing the turbulent magnetic flelds employed 
here. Our main results are presented in §3 for sev- 
eral cases having different properties of the power 
spectrum of magnetic fluctuations, and a comparison 
with measurements is done. Finally, our concluding 
remarks and conclusions are presented in §4. 

2. MODEL 

We describe here the model used to study the 
motion of charged particles {0~^ pickup ions) in a 
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Fig. 2. Configuration used to study the motion of pickup 
ions, of mass m and charge q. The X-axis points toward 
the Sun, the Z-axis points perpendicular to the ecliptic 
and the V-axis completes a right-hand triad. The draped 
magnetic field around Mars goes from a minimum to a 
maximum at the equator and pole, respectively. The 
SW velocity rises from a minimum to a maximum at 
the equator and pole, respectively. The angle 8 is the 
solar zenith angle (SZA). The radius of Mars is taken to 
be Rm = 3397 km. Only the "northern" hemisphere is 
considered in this work. 

stationary background magnetic field with a fluctu- 
ating magnetic and electric field components. The 
magnetic field configuration and the SW flow prop- 
erties are taken to represent, in a first approxima- 
tion, conditions present in the ionosheath of Mars. 
The geometry of our model is shown schematically 
in Figure [21 with several quantities identified below. 

2.1. Equation of Motion 

We follow the three-dimensional motion of par- 
ticles in a time-dependent fluctuating magnetic field 
and its corresponding induced electric fluctuation. 
We will refer to such fluctuations as "turbulence", 
meaning by this that it has a particular power spec- 
trum density distribution as discussed below. We 
treat only the non-relativistic motion of the pickup 
ions. We consider that the dynamics of these ions is 
determined solely by an electromagnetic interaction. 

The equation of motion for a particle of mass m 
and charge 5, in a region with electric E and mag- 
netic B fields is given, in the non-relativistic limit, 
by the Lorentz equation (e.g. Jackson 1975): 



— (mv) = g E + - X B 
at V c 



(1) 



For our problem at hand, the magnetic field is as- 
sumed to be of the form: 



where Bq = Bo(r) is non- homogeneous stationary 
background field, while 6'B{r,t) represents the tur- 
bulent part; with a zero-mean average value. No 
background electric field is considered here, so only 
that E resulting from the fiuctuating 5B will be 
taken into account in the equation of motion (see 
below). 

2.2. Steady Background Fields 

The steady fields determining the dynamics of 
the pickup ions are a non-homogeneous background 
magnetic field Bo and the streaming SW velocity 
Uo ; see Figure [2] 

We are interested mainly in the dynamics of the 
ions over the terminator, hence restrict our simu- 
lation region to a radial interval and y-coordinates 
given by, respectively: 



r e Rout], and y G 



Rp 



R„ 



where = 



Rp is the radius of the planet 



and Rout is an external boundary; to be defined be- 
low. The regions 2; > and z > are the only ones 
considered for our numerical integrations. 

The non-homogeneous background magnetic field 
used here resembles that present in non-magnetic 
planets. The magnetic field Bq in the F-direction is 
taken to be essentially constant, but with a depen- 
dence on the coordinates {x,z). For simplicity, we 
adopt the following analytical expression (RAPIO) 
for this magnetic field, in component form: 



Br 





— B,nin ) cos 9 





B„ 



(3) 



where cos 9 = x/r, and {Bj^in, B^ax} correspond to 
values of the magnetic field a the polar {9 = tt/2) and 
sub-solar points {9 = 0), respectively; see Figure [5J 

The velocity field Uq of the SW plasma around 
the planet is taken to increase sinusoidally in mag- 
nitude from the equator to the pole as: 



— [{Umax — Umin) siu 6* + Umin] siu 6* 



+ [{Umax - U„ii„)sm9 + Ujnin] COS 9 



B{r,t) =Bo + SB{r,t), 



(2) 



(4) 

where {Umin, Umax} corrcspond to values of a sub- 
solar points {9 = 0) and polar {9 = tt/2) points, 
respectively, and sin 9 = z/r. 

2.3. Turbulent Magnetic Field 

There is at the time being no general theory of 
turbulence in plasmas and it appears to be no unique 
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way to describe it (e.g. Borovsky & Funsten 2003, 
Cho et al. 2003, Zhou et al. 2004, Horbury et al. 2005, 
Galtier 2009). On other hand, for example, the value 
of spectral index of the SW turbulence being of a 
Kolmogorov or Iroshnikov-Kraichnan type or its na- 
ture is still a matter of investigation (e.g. Bruno & 
Carbone 2005, Ng et. al. 2010). 

In the case of homogeneous turbulence one may 
construct the magnetic field fiuctuation (5B at a 
particular spatial point by adding the contribution 
of a large set (iVm) of plane waves with differ- 
ent wave-number (k). This can be done, for ex- 
ample, by means of a Fourier transform on P{k) 
(e.g. Owens 1978) or by orienting randomly in space 
the wave-vector number (e.g. Giacalone & Jokippi 
1994), with the field amplitude satisfying a particu- 
lar form of the power spectrum. 

Both approaches above have their pros and cons 
(e.g. Casse et al. 2002) and have been used in dif- 
ferent works on the effects of turbulence. For ex- 
ample, in cosmic ray dynamics in the radio lobes of 
galaxies (e.g. Fraschetti & Melia 2008, O'Sullivan et 
al. 2009) or the acceleration of particles in shocks 
(e.g. Muranushi & Inutsuka 2009). Here we con- 
sider the propagation of one-dimensional turbulent 
plane waves traveling in each perpendicular direction 
in turn, and assume homogeneous turbulence. This 
will allow us to assess the effect of maximum tur- 
bulence effects along each direction on the dynamics 
of the ions; which is the main objective of this 
work. 

We write the magnetic fiuctuation field at a par- 
ticular point in space as a superposition of trans- 
verse plane waves propagating, for example, in the 
±Z-direction. Hence we write a Z-propagating 
wave with arbitrary polarization (e.g. Jackson 1975, 
Damask 2005) as the real part of: 



The real value functions 6b j = 6b{kj) in ([5]) are 
the fluctuating field amplitude corresponding to the 
j-th mode of the m-number of modes Nm- The fre- 
quency of oscillation of each mode is ujj = ivphkj, 
with Vph being the phase-velocity of the j-th mode. 
For definitivencss, the phase-velocity of the pertur- 
bations is taken here to be the Alfven velocity of 
the plasma (i.e., Vph = va)- The sign in the ex- 
pression for LJj accounts for waves traveling both 
in the positive or negative Z-direction, respectively. 
The form of turbulence represented in ([5]) satisfies, 
by construction, the divergence-free character of the 
magnetic field fluctuation (V • (5B = 0). A similar 
construction is done for each orthogonal propagation 
direction considered here. 

The amplitudes Sbj in ^ are chosen to satisfy a 
particular one-dimensional power spectrum of fluc- 
tuations P{k). We assume here that the power spec- 
trum density P{k) is given by a power-law: 



P{k) = Pjv {SB^) k-^ 



(7) 



where {SB"^) is the magnetic held fiuctuation vari- 
ance, 7 the spectral index and Pn a normalizing con- 
stant. For 7 = 5/3 we have a Kolmogorov (1941) tur- 
bulence and for 7 = 3/2 the Iroshnikov-Kraichnan 
type spectrum (Iroshnikov 1963, Kraichnan 1965). 
The normalization of the spectrum is such that 



/•OO 

/ P{k) dfc = {5B^ 
Jo 



(8) 



In practice, the limits of integration in ([8|) go from 
a kmin to a kmax dictated essentially by numerical 
criteria of the calculations. 

The discrete amplitudes 6bj in ([5]) can now be 
obtained from the power spectrum ([7]). Using the 
orthogonality of the polarization state vector, 



5B 



6B, 
SB, 



Y,Sb,\p)e' 



(kj z—cujt) 



(5) 



where the polarization state or Jones vector \p) in 
the plane is given by 



\P) 



Px 




cos If e"^=^ 






sin (f c'*!' 



(6) 



with the (/s's and </)'s being parameters describing 
the polarization state of the wave. For example, for 
02; = 01/ we have a linearly polarized wave inclined 
an angle ip, while if (jj^ and 4>y differ by 7r/2 circular 
polarization is produced; otherwise elliptical polar- 
ization is achieved. 



(p\p) = [pI^pI] 



Px 

Py 



where {p\ is the complex conjugate (dual) vector of 
\p), and {SB"^) ~ SB ■ SB^ and assuming statistical 
independence among the different wave modes, we 
have 



N„ 



Y,p{k,)\k, 



Sb) = P{kj)Akj 



(9) 

The numerical realization of our turbulent mag- 
netic field is obtained as follows. Given a mean fluc- 
tuation ((5_B^), a spectral index 7 and {kmin,kmax} 
for the fluctuations, a normalization constant Pn is 
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Fig. 3. Intensity plot of a particular turbulent magnetic field realization of plane waves traveling in the Z-direction at 
different times. The X-component is shown at the left figure while the Y-component at the right figure. The total field 
variance here is {5B^) = f and 7 = 2. Units are arbitrary. 



determined from ^ . We choose a set of Nm equally 
spaced wave-numbers in log /c-space in the relevant 
interval. From this set of k values, the amplitudes Sbj 
are obtained using ©. Since we are not interested 
in any particular effect of the polarization state of 
the turbulent waves, we choose random phase values 

and (jjz and angle in ((5]) uniformly distributed 
in [0, 2tt] for each kj value. 

At each point, for example, z and time t the sum 
in ([5]) is evaluated, and the SBx and 6Bz components 
are obtained. In this form, the total magnetic fluc- 
tuation (5B is constructed; which is to be added to 
the inhomogeneous background field Bq of equation 
^ . A similar construction is done for each propaga- 
tion direction of turbulent waves when required. In 
Figure [3] we show a particular random realization of 
a trave;omg turbulent field. 

2.4. Electric Field 

In the stationary reference frame the electric field 
depends on the local total velocity field of the SW 
plasma U around the planet and the local magnetic 
field. In the non-relativistic limit this leads to 



E(r,t) 



1 



U X B(r,i) , 



(10) 



where B is the total magnetic field obtained from 
equations ^ and ([3]). The total velocity field U is 
obtained as U = Uq + SU, were we consider that 
dXJ{r,t) ~ —vaSH/Bq as in a MHD Alfvcnic wave. 

The fields determined by equations ([2]) and ([T0|) 
are used in equation ([l} to determine the dynamics 
of each pickup ion in the ionosheath of Mars. 

2.5. Initial Conditions and Integrator 

For a particular mode in the plane-wave expan- 
sion ([5]) one expects that if the time-scale 



for the field fluctuation to be about the same as the 
time-scale of gyrotropic motion of the ions (~ 
with n ~ qB/mc), the dynamics of the particles will 
drastically be affected. On other hand, in order for 
an Alfvenic-type turbulence to affect the motion of 
pickup ions it has to modify such motion in the time 
scale spent for the ion to transverse our region of 
interest; e.g., from the subsolar point to the termi- 
nator. These considerations guide us, as well as in 
situ measurements, to adopt some numerical values 
for different parameters in our calculations. 

As our fiducial model for the ionosheath of Mars 
we adopt the following set of values for different 
parameters in our model. The range in magnetic 
field is taken as {Brmn,Bmax) = (4, 7)Bo, where 
Bq = 5 nT, that we choose based on the data from 
Bertucci et al. (2003). The SW plasma velocity is 
taken as Vsw = 400km/s and the Alfven velocity 
of waves va = O-^Vsw- The range of velocities of 
the SW plasma flow around the planet is assumed to 
be {Urmn.Umax) = (0.01 , 0.5) , consistcnt with 
values measured in that region; see for example the 
data in Figure [TJ Under these conditions, the mean 
local cyclotron frequency for a proton in our region 
of interest is /p = Q.p/2TT = 0.42 Hz; where the mean 
field is taken to be {B) = 5.5i?o- The range of fre- 
quencies for the turbulence spectrum is taken from 

frmn = /p/16 « 0.03 Hz tO fmax = ^Ofp ~ 13 Hz, 

that is consistent with the range in P{k) measured 
in the ionosheath of Mars (e.g., see Figure 5b in 
Grebowsky et al. 2004). For the va chosen, the 
following range in length-scales (A = 27r/fc) follows 
iXrnin,Kiax) « (3.2, 1526) km. 

The turbulence spectrum index is taken to be 
7 = 2 in concordance with measurements in the 
ionosheath of Mars. The amplitude of the turbulent 
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component of the magnetic field is parameterized by 
the energy ratio e= {SB^)/Bq, and we have consid- 
ered the particular values e = (0.0, 0.1, 0.5, 1.0). 

The equation of motion ([T]) for each ion was 
solved using a fourth order Runge-Kutta algorithm 
with an adaptive time-step instead of the standard 
Buneman-Boris algorithm (e.g. Birdsall & Lang- 
don 1985), a time-symmetric second-order scheme, 
used in plasma dynamics. This was adopted after 
several numerical tests with both integrators, and 
by the requirement of following more accurately the 
trajectories of particles especially when a large am- 
plitude of turbulence (e w 1) and rapid oscillations 
were present in the Martian ionosheath. 

A total of 10^ charged particles were used for 
each value of e and propagation direction of turbu- 
lent waves considered, yielding a total of one million 
particle integrations in our numerical experiments. 
The 10^ particles were obtained from sets of 10"^ 
particles evolved in 10 different random realizations 
of the power spectrum P(fc), using each time a dif- 
ferent initial seed for the random number generator 
required to construct the turbulent 5B. Particles 
were distributed spatially only in the "upper" day- 
side of Mars, and following an exponential density 
decay with altitude as in RAPIO. 

3. RESULTS 

In this section we present results obtained for the 
motion of ions under different conditions of tur- 
bulence. We first present the space trajectories of 
the ions and afterwards their velocity profiles at the 
terminator. 

3.1. Ion Trajectories 

In Figure U we show twenty random trajectories 
of O"*" ions in the XZ plane for turbulent waves prop- 
agating in the x, y and z-direction for the extreme 
amplitude e = 1, along with the e = case, and in 
Figure [5] the corresponding ones in the YZ-p\ane. 
The square box in these figures correspond approxi- 
mately to an altitude range from 300 km to 2300 km; 
see Figure [TJ 

The motion of pickup oxygen ions without 
turbulence (e = 0.0), is qualitatively similar to 
that found in previous works on the dynamics of 
pickup ions with an analytical model for the B- 
field (e.g. Luhmann & Schwingcnschuh 1990, Luh- 
mann 1990, Kallio & Koskincn 1999, RAPIO) and 
from that derived from MHD simulations (e.g. Jin 
et al. 2001, Bofiwetter et al. 2007, Fang ct al. 2008, 
Fang et al. 2010). 

The presence of even a small amplitude (e = 0.1) 
of turbulent waves propagating in the ionosheath 
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Fig. 4. Trajectories of O"*" ions in a model of the 
ionosheath of Mars in the XZ-plane. Only a turbulence 
power of e = and 1 has been considered here. The di- 
rection of propagation of turbulence of the wave-vector 
k is indicated at each panel. 
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Fig. 5. As in Figure [4] but for motion projected on the 
Y Z-p\&ne. 

shows a strong effect in the yZ-motion of pickup 
ions, while it is less noticeable in the ATZ-plane. 
Even in the case of e = 1.0 the XZ-motion is not 
greatly affected (FigureS]), in particular the trajecto- 
ries are not importantly "bent" toward the SW direc- 
tion of motion as is inferred from the Mars-Express 
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Fig. 6. Phase-space diagrams of O"*" ions at the Martian terminator: Vz — z and \Vx 
in the X-direction of increasing amplitudes e are shown from left to right. 
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measurements (Figure [T]). 

The dynamics of pickup ions is clearly very com- 
plex even under a small "perturbation" to the back- 
ground field that, for example, may preclude any es- 
timation of the ionosheath conditions based solely in 
the analysis of the motion of pickup ions. The anal- 
ysis of the statistical motion of the ions is out of the 
scope of this paper, but we expected that they will 
have a dynamics similar to that of a Levy random 
walk process (e.g. Shlcsinger, West & Klafter 1987, 
Greco et al. 2003). 

3.2. Velocity Profiles 

Several processes may be investigated from the 
velocities of the ions in the turbulent medium, such 
as their phase-space structure or the pitch-angle 
difussion and scattering (e.g. Price & Wu 1987, 
Li et al. 1997, Cravens et al. 2002, Fang et al. 2008). 
Although important for our understanding of 
pickup ions, our focus here is essentially in the ver- 
tical velocity profiles of the ions over the terminator 
in order to compare with the data shown in Figure [T] 

The vertical velocity profiles of ions were 
computed at the terminator, and considering only 
a small region of width in the y-direction of Ay = 
±i?M/4; as indicated in Figure O Twenty bins in 
the Z-direction were set and the average velocity of 
all particles crossing the terminator were computed 
for each component. This was done for each level of 
turbulence and direction of the Alfvenic waves. 

In Figure H] we show "phase-space" plots of the 
particles above the Martian terminator, in the box 
described above, for our fiducial parameters. The 
different graphs correspond different levels of turbu- 
lence (e = 0,0.1, 1) propagating in the X-direction. 
A clear effect of heating and scattering in phase- 
space due to turbulence is appreciated as its am- 
plitude £ increases in value. Under our fiducial 



values, and the different prescriptions for turbulent 
waves, we notice already from Figure IH] that the Z- 
component of the velocity dominates over the X- 
component. 

In Figure [7] we show the average velocities of Vx , 
Vy and Vz of the ions, for turbulence amplitudes 
e = 0.1 {top) and e = 1.0 {bottom). The case of 
e = has a similar behaviour as that of e = and 
we do not consider it in Figure [T) The graphical 
representation is similar to that in Figure [1] in order 
to aid in the comparison with the in situ data; the 
Z-component of the velocity has been inverted in 
sign for that purpose. 

The trend displayed of the velocity profiles is 
similar to that reported in RAPIO, although signif- 
icant differences appear now when turbulent waves 
are taken into account. In the no-turbulent case, 
e = 0, the (vy) is zero since no motion along the 
y-direction is generated (see Figure [5]). Once a tur- 
bulent component sets in, the erratic motion of par- 
ticles contribute to a non-zero value of (vy). 

A clear result from all of our numerical calcula- 
tions is that even under a strong (e = 1) level of 
turbulent waves propagating in any direction, the 
trend of the Z-component of the velocity of ions al- 
ways dominating over the X-component is preserved, 
contrary to measurements at such altitudes in the 
Martian terminator. 

Unfortunately, for example, the extensive and 
complete MHD simulations of Fang et al. (2008) do 
not address the velocity profiles of pickup ions at the 
terminator, so we are not able to compare our results 
with their work. In this respect, it will be interesting 
to measure from this kind of MHD simulations the 
vertical profiles of 0+ ions and compare them with 
data. 
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Fig. 7. Vertical velocity profiles of O ions at the Martian terminator. From left to right, the plots correspond to 
different turbulent waves propagation vector. The upper row correspond to values of e = 0.1 and the lower for e — 1.0. 



4. DISCUSSION AND FINAL COMMENTS 

4.1. Discussion 

We have calculated the trajectory and average ve- 
locity profiles of newly born ions in the dayside 
ionosheath and magnetic polar regions of Mars, as 
they are accelerated by the convcctive electric field 
due to the streaming solar wind plasma. In addi- 
tion to a simplified model of the draped IMF, which 
determines the large-scale geometry of the magnetic 
field, we have included for the first time an approx- 
imation to the "turbulent" component for the mag- 
netic field in that region, and studied its effect on 
the dynamics of pickup ions. We have compared the 
velocity profiles resulting from our calculations to 
those measured in the region with the ASPERA-3 
instrument onboard the Mars Express, as reported 
by PdT09. 

The effect of a turbulent part (5B in the magnetic 
field is most significant in deflecting the trajectories 
in the plane perpendicular to the SW flow (plane 
YZ), but it has little effect on modifying the ve- 
locity structure of pickup ions at the terminator in 
comparison to the absence of turbulent waves. This 
behavior holds under different power-laws (7) for the 



turbulence power spectrum and its amplitude (e), as 
well as for its direction of propagation and small vari- 
ations of our fiducial values ( ij2.5p of the parameters 
in the simulations. 

In searching for a pure dynamical explanation of 
the velocity profiles measured by the Mars-Express 
spacecraft, we decided to use parameters in our 
model that would probably favor the dominance of 
the Vx component over the Vz component. For 
this, we would require a strong magnetic field at 
the poles, Bmin, and a small velocity Umax, since 
the gyro-radius of a particle, for the same q/m ra- 
tio, is rg oc v±/B and this would make the trajec- 
tories bend more such that at the terminator the 
X-velocity would be increased. 

We decided to use for a new set of runs with 
Bmin = GBq = 30 nT and Umax = O.S75Vsw = 
150 km/s, and keeping the rest of the parameters the 
same as in our fiducial case. As one might expect 
there is a dominance of the Vx component over the 
Vz , but this occurs only for certain altitudes as shown 
in Figure |S] and [HI Even with these values of the pa- 
rameters, rather extreme according to data (see, for 
example. Figure [1] from where an estimate of Umax 
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Fig. 8. The same as in Figure [6] but with parameters more favorable for a strong curvature in the paths of O"*" ions. An 
increase of particle density in different regions of the phase-space is noted in comparison to that ofFigure [B] 
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Fig. 9. As in Figure [7| but with a stronger Bmin = 30 nT and lower Umax = 150 km/s in order to favour an increase 
in the curvature of the pickup ions and its Vx component. Note that the general behaviour of the profiles do not 
correspond to the ASPERA-3 measurements. 



can be obtained) the behaviour of the profiles do not 
correspond to the measurements; for example, no in- 
version of the velocity profiles at altitudes ~ 1000 km 
is present. This numerical experiment strongly sug- 
gests that even under more propense conditions, al- 
beit of not being consistent with measurements of 
the properties of the i3-field and stream velocity of 
the SW at the pole, the idea of the behaviour of 



being due to charge-particle dynamics is not sustain- 
able. 

On the basis of our results, we conclude that 
turbulent-like magnetic field fluctuations, and the 
accompanying convective electric field variations, are 
not likely to influence the dynamics of picked-up 
ions in a sufficient manner to explain their essen- 
tially tailward velocity measured over the magnetic 
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pole of Mars. Furthermore, the relative behaviour 
of the different components of the velocity obtained 
in particle simulations does not correspond to what 
is measured at the Martian terminator. This sug- 
gests that additional processes to a purely charged- 
particlc dynamics must be playing a dominant role 
in the interaction of the collisionless solar wind and 
the ionospheric plasma in the region. 

A possibility is that advocated by Perez- 
de-Tejada and collaborators (e.g. Perez-de-Tejada 
& Dryer 1976, Perez-de-Tejada 1999, PdT09, 
RPAVIO), where a viscous- like interaction of the 
SW plasma and the ionoshcath plasma is occur- 
ring. The interaction of the SW with the upper 
ionosphere over the magnetic polar regions of non- 
magnetic planets and objects is characterized by the 
existence of several features that lend themselves to 
be explained in terms of a fluid viscous-like interac- 
tion. Several of the phenomena for which it provides 
an explanation are briefly summarized here for the 
sake of completeness. (1) The existence of a veloc- 
ity shear in the ionosheath as spacecraft approach 
Venus over the magnetic poles (Romanov et al. 1979, 
Bridge ct al. 1967, PdT09). The accompanying de- 
crease in plasma density and the increase in the 
gas temperature are consistent with a viscous-like 
boundary layer in the region. (2) A similar situa- 
tion is observed in the ionosheath of comet Halley 
below the so-called "mystery transition" , where the 
decrease in SW velocity and density and the simul- 
taneous increase in the gas temperature, measured 
by the Giotto spacecraft, is consistent with a viscous- 
like flow characterized by an effective Reynolds num- 
ber «30 [Perez-de-Tejada 1989, RPAVIO]. (3) The 
transterminator flow observed by the Pioneer Venus 
Orbiter (PVO) is consistent with a viscous-like drag- 
ging of the Venus upper ionosphere by the SW flow 
(Perez-de-Tejada 1989). Since such interaction oc- 
curs preferentially over the magnetic polar regions 
this leads to the carving-out of plasma channels in 
the ionosphere and near wake, providing a simple 
explanation for the ionospheric holes measured by 
PVO in the nightside of the planet (Perez-de-Tejada 
2001 & 2004). It is important to indicate that the 
physical nature of such "anomalous viscosity" is not 
addressed in these papers, and is the subject of cur- 
rent investigations. 

4.2. Final Remarks 

It must be stressed that the calculations pre- 
sented in this paper are not fully self-consistent, in 
the sense that the electromagnetic field and the so- 
lar wind velocity around the Martian ionosheath are 



assumed, not solved for. Nonetheless, the assumed 
SW velocity field and i3-field structure are represen- 
tative of what might be the Martian magnetic field 
in the region of interest; the case of the southern 
hemisphere would be more complex to model due to 
the presence of inhomogenictics due to crustal fields. 

The type of approach taken here clearly prevents 
us from studying plasma instabilities which may play 
an important role in the coupling of solar wind and 
ionospheric plasma and significantly change our con- 
clusions. The effect of such instabilities on the dy- 
namics of ions, as well as that of magnetic field 
features related to its 3D character and the pres- 
ence of magnetic concentration regions, are beyond 
the scope of the present study and will be treated in 
future contributions. 

This research was funded by UNAM-PAPIIT 
Research Projects IN121406 and IN109710, and 
CONACyT Projects 25030 and 60354. 
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